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Abstract 

We study a simple model of the two-channel Kondo effect in a distorted 
crystal. This model is then used to investigate the interplay of the Kondo 
and Jahn- Teller effects, and also the Kondo effect in an impure crystal. 
We find that the Jahn- Teller interaction modifies the characteristic energy 
scale of the system below which non-Fermi-liquid properties of the model 
become apparent. The modified energy scale tends to zero as the limit 
of a purely static Jahn- Teller effect is approached. We find also that the 
non-Fermi-liquid properties of the quadrupolar Kondo effect are not stable 
against crystal distortion caused by impurities. 



1 Introduction 

The general multichannel Kondo model describes a magnetic impurity of spin 
S in a sea of conduction electrons having spin-^ and n degenerate orbital 
channels^. Such a model was found by CoxQ also to apply to a nonmagnetic 
rare earth impurity such as U 4+ or Ce 3+ with an orbitally degenerate doublet 
ground state. In this case, a pseudospin variable associated with quadrupolar 
deformations of the local orbital plays the part of the impurity spin in the orig- 
inal multichannel Kondo model, and the real magnetic spin of the conduction 
electrons plays the part of the channel index. Thus we study the n = 2 channel 
Kondo model. This is known as the quadrupolar Kondo effect. 

This realisation of the multichannel Kondo model is particularly interesting, 
as the spin symmetry of conduction electrons ensures channel symmetry in the 
model. The channel-symmetric two channel Kondo model has non-Fermi liquid 
low temperature behaviour, shown by logarithmic divergence of the impurity 
magnetic susceptibility and specific heat coefficient [0, ||[ 

In this paper we study the interplay of the quadrupolar Kondo effect with 
local crystal distortion. We consider two possible sources of this distortion: the 
Jahn- Teller effect^] and local distortion caused by non-magnetic disorder in 
the crystal. For an ion subject to the quadrupolar Kondo effect, the orbital 
degeneracy of the ion realises conditions for the Jahn- Teller effect: an orbital 
degeneracy will split spontaneously on coupling to phonons. This splitting is 
accompanied by a distortion of the crystal around the ion. In this case, there 
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will be several distortions producing the same splitting of the orbital degeneracy, 
and the system may tunnel between these equivalent distortions (the dynamic 
Jahn- Teller effect). There has been discussion[Q, ||, [7| of this interplay of the 
Jahn- Teller and Kondo effects. It is the intention here to provide a simple model 
of this interplay, and use the technique of Abelian bosonization to calculate the 
low-temperature properties of this model. 

Also, other (non-magnetic) impurities in the crystal will cause internal 
strains. The effect of strain on a quadrupolar Kondo ion may be described 
by a term like Hi = gs z in the Hamiltonian, where g is linearly related to com- 
ponents of the local strain tensor: i.e., it enters the Hamiltonian like a magnetic 
field term. We investigate the effect of internal strains on the non-Fermi liq- 
uid behavior of a quadrupolar Kondo impurity, since this behavior is the main 
motivation for studying such systems. 



2 The model 

Cox|| has shown that many different crystal symmetries are sufficient to pro- 
duce a quadrupolar Kondo effect in impurity ions (see also Hirst || for a review 
of coupling between a local moment and conduction electrons). Since we are 
studying the interplay between this effect and local distortions, we choose a 
simple crystal symmetry with the hope of producing a tractable model. Thus 
the model we consider is one of an impurity ion at a site of tetragonal sym- 
metry (lattice lengths a = b 7^ c). The impurity is taken to be subject to a 
quadrupolar Kondo effect. The three-dimensional Kondo Hamiltonian may be 
reduced to a one-dimensional Hamiltonian by partial- wave projection (see Ref. 
II Appendix A for a review of this), which in the two channel case is 



(1) 



H K = J2 [H m + ^+J~(0) + H.c.) + l8 t J?(0)} > 
where 

J-(x) =■ V^(x)^ T (x) : and Jf(x) = - : xp^ (x)^ (x) - ^(aj)^;(x) : (2) 
are the spin currents in the two channels and 

/+00 
dx : 4 jtJ (x)(-id x )^ a (x) : (3) 
-00 

c=ia 

describes right-moving fermions on —00 < x < +00. This is discussed in the 
review by Cox and Zawadowski(l(J. 

To this is added a part describing local distortion. The Jahn- Teller distor- 



tion may be described jTTJ by a pseudospin r. In the current case, only the 
distortion of the four nearest neighbours in the a6-plane are considered, as only 
this distortion will couple to the orbitals giving rise to the Kondo effect. Thus 
we may take t z = \ to correspond to the distorted state with an expansion 
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in the x-direction and contraction in the y-direction. Likewise t z = — ^ cor- 
responds to the state with expansion in the y-direction and contraction in the 
x-direction. 

These distortions are coupled to the pseudospin s in via a term 

Hi = gs z r z (4) 

describing the splitting of the degenerate orbitals by the Jahn- Teller distortion. 
A further term 

H 2 = A t x (5) 

is added to describe tunnelling between the states t z = ±i. The total Hamil- 
tonian is then 

H = H K + Hi + H 2 . (6) 

Note that to describe the effect of local strain caused by non-Kondo impurities, 
we set Ao = 0. g is then linearly related to components of the local strain 
tensor. 

We now follow Emery and Kivelson[|j and use Abelian bosonization and 
refermionization in the analogue of the Tolouse limit for the one-channel Kondo 
model to transform our Hamiltonian ^ to a more tractable form. Bosonizing 
via 

^ a (x) = ^L^v^M (7) 
and introducing the spin and charge combinations of the fields 

<Pci(x) = -j=[(j)^{x) +<t>il(x)] 

<j>i(x) = —[faffl-fa^x)] (8) 
and the spin and spin-flavour combinations of the spin fields 

<t>s( x ) = ^g[0l( a O + #2( a O] 

0sfO) = -j=[<j>i{x) - fc{x)] (9) 

we obtain 

H K = H [(j> s ] + H o [0 s f] + -^( a+ ^M0) cos[v^0 s f (0)] + H.c.) 

+ -L Sz d x <p s (o) (io) 

'IT 



where we have dropped, the charge fields <p c i 

as they decouple from the impurity. 

Using the canonical transformation 
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H K ^U^H K U = iT o [0s] + #o[0sf] + — s, cosfV^^f (0)] 

7ra 

+ ^± Sz d x cP s (0) (12) 



with A+ = I — 2~KVf. 

We take A+ = in the following, analagous to taking the Tolouse limit in 
the single-channel Kondo problem. Deviations from A+ = are discussed in 
section 3, and are irrelevant. This has the effect of decoupling the <p s field from 



the impurity. Defining for notational simplicity <j)(x) = v47T0 s f (a;) and A = — 
we have the bosonised Hamiltonian 

H = H [<j)} + Xs x cos[<£(0)] + gs z r z + A t x (13) 

The next step is to refermionise this Hamiltonian: this is similar to the 
refermionization procedure used by Moustakas and Fisher [12] in their study 
of two-level systems in metals. Care must be taken with the fermionic repre- 
sentation of the pseudospins s, r in order that the refermionized Hamiltonian 
be quadratic and that all commutation relations are satisfied. The standard 
fermionic representation of spins s*i , s 2 is 



sx+ = d\ , si z = d[di - \ 

S2+ = 4j*4*i , S2z = d\d 2 - \ 



(14) 

and these components are rotated before identifying s i <-> f and s 2 «-> s so that 

(15) 






Then refermionizing </> via 

e 1 ^ = V2^e in ^ dl+d ^ d2 ^(x), (16) 
gives the quadratic Hamiltonian 

H = H [il>] + ^V2^(4-d 2 )[^(0) + i;(0)} 

+ i(di-d 2 )(d\ + d 1 ) 

+ Ao(4di-±) (17) 

where Hq[iP] describes free right-moving fermions on— oo < x < +oo. 

The factor exp[i7r(d{cZi +d\d 2 )] introduced in Equation (|Ul) serves to ensure 
the commutation relations of the spins with the electron fields are satisfied. It is 



precisely the Klein factor discussed in the review by von Delft and Schoeller[13 
Transforming to Major ana fermions via 

i/i(x) = —r=[x(x) +i£(x)\ 
d i = "^f^ + ic) 

d 2 = ^W + ifi) (18) 
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the £-field decouples from the impurity leaving 

H = H [x] + iAax(O) + igab + iA 6c 
l - -~ A A — Ac 



where A = — y ^A, g = — ^g and Ao = /ao 
Impurity Green's functions, defined as 



(19) 



Ai(t) = -i(TCi(t)Cj(0)> 



(20) 



where £ = (o,,b,c), can be calculated from the Hamiltonian (19) by e.g. the 
equation-of-motion method. In Matsubara formulation we obtain 



Dij(iu n ) 



1 



A 2 G(°)( 



(21) 



\UO r . 



A 2 



-gun 

3A0 



-iw„[io;„ - A 2 G(°)(io; n )] -iA [iw„ - A 2 G(°)(kd n )] 
iA [iw„ - A 2 G(°)(i W „)] co 2 + 5 2 + k, n A 2 G(°)(iw n ) 



where 



G (o)( 



2w F 



signa; r . 



(22) 



is the free conduction electron Green's function. 

The impurity contribution to the free energy F may be calculated by evalu- 
ating the thermodynamic average of ax(0) using ( |2l| ) and integrating over the 
coupling constant: 



5F 



iA / da(ax(0)) 
Jo 

~TJ2 I'd* 



Tsigna; n 



W n , .2 1 a ^ r CKl Sign (J n 



(23) 



where T 



At 

2v F 



The sum over uj n may be changed to an integral in a standard manner: 



F = F Q + 



duj 
2^ 



f(u>) tan 1 



(24) 



where -Fo is the free energy of the impurity decoupled from conduction electrons. 
It is necessary to keep a momentum cutoff $7 to obtain a finite expression 
for F, but may be taken to infinity in calculating quantities obtained by 
differentiating F. 



3 The Jahn- Teller distortion 

The properties of model ( |l9| ) in the case Ao 7^ (modelling Jahn- Teller distor- 
tions) are very different from those in the case Ao = (modelling deformations 
caused by local strain fields). This section is thus devoted to the case Aq 7^ 0. 
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From (|24j) the impurity contribution to the entropy may be obtained 

c _ 9F 
b - dT 

— tttt tan 



= 5 -<5S (25) 

where So is the entropy of the impurity decoupled from conduction electrons. 
This can easily be calculated as 

1 W &- E +/ T _i_ P „-JS_/T 

S„ = ln(2) 4- N.-^ 4- + l gt^_±|^L_ (26) 



where -E-j- = ±J g 2 + An, which has a low temperature expansion 



So = ln(2) + + l) e~^ T + O (^e"**^) . (27) 

(55 1 may be expanded at low temperature by using 
/ 



tan 1 



vr -, / Too g 2 + A2\ /, ,o 

r , - ^ sign W - tan" 1 !— ) + Q f (r W ) 3 



(28) 



to obtain 



2 V ; 6r A 2 



g = lin( 2 ) + ^- a +0(T 2 ). (29) 



Equation ( |2 X|) may also be used to examine the local impurity susceptibility 

X i(T)= ^ dr((s z (r)s z (0))), (30) 



o 



i.e. the susceptibility to a strain field acting on the impurity only. From the 
refermionization prescription in Equations (|14D and (|T8D, we have 



s z = iV2abc (31) 

and thus ((s z (t)s z (0))) must be expanded using Wick's theorem: 

((8*(t)*z(0)» =-2 [^(O)^^)^^) - ^6(0)^^)^(0) 

+Z> a6 (0)L> cc (r)A l6 (0) - D ac (0)D ba (r)D bc (0) 
+D ac {0)D bb {r)D ac {0) - D ac (0)D bc (r)D ab (0) 
+D aa (T)D bc (0)D bc (0) - D aa (T)D bb (T)D cc (r) 
+D aa (T)D bc (T)D cb (T) - D ab (T)D bc (0)D ac (0) 
+D ab (r)D ba (T)D cc (T) - D ab (r)D bc (T)D ca (T) 
+D ac (r)D bc (0)D ab (0) - D ac {T)D ba {r)D cb (T) 
+D ac (T)D bb (r)D ca (r)}. (32) 
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Since we are interested in the low-temperature properties of the model, 
we find the major contributions to this by finding the discontinuous Green's 
functions, Equation (21). At zero temperature, a discontinuous Dij(\uj n ) will 
give a Dij{r) behaving as 1/r at large r; whereas a continuous Dij{iu n ) will 
give a Dij (r) behaving at worst as 1/t 2 . We find that D aa , D ac = D ca and D cc 
are discontinuous in the ^-representation. Then defining the constants 

Ai = D ab (0) 

A 2 = D ac (0) (33) 
^3 = D bc (0) 

we find 

((**0X(0))) ~ -2[A 2 3 D aa (r) + 2>M3A»(t) + A 2 D cc (t)], (34) 



these terms arising from terms 1, 3, 7 and 13 of Equation (p2\j. 

We obtain large-r properties of this by examining small-u; properties of the 
Dij(\u n ). Since 

DaailUJn) = 7^ sign LJ + 0(u) 

D ac (iuj n ) = — ^— signw + 0(u) (35) 
ll A 

g 2 

DccQwn) = — — 2 signw + O(w) 
11 A Q 



we have at low temperature and large r 



1 



7rrr 



M^n) ~ (36) 



-Dcc(iw ri 



5 2 



ttTA 2 t 



Thus 



2 / g „ \ 2 /-l/T-ro dr 



r 



^ + zHM4> (37) 



where tq is a cutoff chosen so that the asymptotics ( pq ) are applicable (i.e. this 
result holds for T <C to). 

We may also study deviation from the Tolouse limit; i.e. A+ 7^ in Equation 
(|l^). The operator we are including is 

o+ = ±± Sz d x uo) (38) 
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and if we refermionize 4> s and transform to Majorana fermions via 



l 



\/2vrao 
1 



'W iircl> s (x) 



we find the following expression for this operator: 

0+(r) = -a(r)6(r)c(r)xs(0,r)| s (0,r). 



(39) 



(40) 



At the Emery-Kivelson line the impurity spin s z is decoupled from the ip s elec- 
trons so the 0+— 0+ correlation function factorises into the spin-spin correlation 
and the density-density correlation function for free electrons. 



«O + (r)O+(0)» = <Mr)*,(0)» [G(°)(7 



(41) 



Since 



1/r 
1/r 



G(°)(t) 
(('.(r)*,(0))) 
at large r and zero temperature, we have 

«O + (t)O + (0)» ~ 1/r 3 



(42) 



(43) 



and thus the operator + has scaling dimension 3/2 and is irrelevant. 

However, as found by Sengupta and Georges [15] it is this operator which is 
responsible for the behaviour of the impurity specific heat. Comparison with 
the calculation in Ref. |l4j gives the result (again for T <C to) 



dC imp (T) 



dT 



2ttFvI 



In 



1 

T7o 



(44) 



Comparing Equation (29) with the T — > entropy of the two-channel Kondo 
model Hi 



S(0) = ± In 2 



(45) 



shows that the Jahn- Teller interaction does not change the residual entropy of 
the two-channel Kondo model. This residual entropy is an important and pe- 
culiar feature of non-Fermi-liquid two-channel Kondo models. Since it persists 
in the presence of a Jahn- Teller interaction, it is clear that the non-Fermi-liquid 
behaviour also persists. 



The results [ 14 [ for the two-channel Kondo local impurity susceptibility 



xf(T) 

and specific heat coefficient 

dT 



In 



A 2 + 



8vrrt;| 



In 



(46) 



(47) 



8 



show that the log-divergence of these quantities is governed by the characteristic 



energy scale of the system, V. Comparison with Equations (37) and (0) show 
that this divergence persists but with a modification of the characteristic energy 
scale to Tq 1 . We find as Ao — > 0, 

n 2 4- A 2 

r ° ~ "TaT (48) 

and thus the non- Fermi liquid behavior emerges only at lower and lower temper- 
atures. This is because the dynamic nature of the Jahn- Teller effect preserves 
symmetry of the system, whereas a purely static Jahn- Teller effect - identical to 
local distortion of the crystal - breaks the symmetry of the system and destroys 
the non-Fermi liquid behavior. 

Thus the Jahn- Teller effect does not destroy the quadrupolar Kondo effect, 
but modifies the characteristic energy scale of the system. 



4 Crystal distortion by impurities 

Setting Ao = in (|^) leaves us with a model for a quadrupolar Kondo impurity 
sitting in a strain field caused by other impurities in the metal. We wish to 
obtain properties of this model which have been averaged over all possible 
distributions of these impurities. Thus the local impurity susceptibility will be 
given by 

(Xi(T)) g = [ P[g] X i(T,g)dg (49) 



where P[g] is the probability distribution of the local strain field, and X i(T,g) 
the impurity susceptibility of the two-channel Kondo model in a field, which 



may be calculated from (24). 

The quantity P[g] is calculated in Refs (|l(|]) and (Jl?]]). It is Lorentzian at 
small g and small impurity concentration, with Gaussian tails at high g. For 
calculational simplicity, we use 

P \9] = -^— 2 i (50) 
7r g z + or 

i.e. a Lorentzian distribution for all g. This is possible as impurities in large 
strain fields have very low susceptibility, so the large-g behaviour of P[g] is 
unimportant. 

Xi{T-,g) has the properties that at g = 0, it is logarithmically divergent as 
T — > 0, but this logarithmic divergence is rounded off at finite g, and Xi(T,g ^ 
0) — > constant as T — > 0. Thus we may use 

X i(T,g) =\n{T 2 + g 2 ) (51) 

to obtain qualitative results on the behavior of (xi(T)) g . 
These approximations give the result 

( Xl (T)) g ^-HH+T) (52) 
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i.e. the impurity susceptibility is finite as T — > 0. Thus the non-Fermi liq- 
uid properties of the quadrupolar Kondo effect are not stable against crystal 
distortion caused by impurities. 

The effect of impurities on the low-temperature behaviour of the impurity 
susceptibility was discussed by Dobrosavljevic et. al. ||l8fl . They considered 
the effect of local density-of-states fluctuations caused by the disorder on the 
susceptibility They found that these fluctuations enhanced the non-Fermi- 
liquid properties of the model. 

The effect we considered above suppresses the non-Fermi-liquid behaviour, 
and is thus in competition with that of Dobrosavljevic et. al. 
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